1 



(N 
> 

o 



> 

m 
o 



(N 



X 



Gauge Invariant Smearing and the Extraction of Excited State Masses 
Using Wilson Fermions at = 6.2 

UKQCD Collaboration - presented by Sara Collins 

Department of Physics, The University of Edinburgh, Edinburgh EH9 3JZ, Scotland 

We present an investigation of gauge invariant smearing for Wilson fermions on a 24^ x 48 lattice at (3 — 6.2. 
We demonstrate a smearing algorithm that allows a substantial improvement in the determination of the baryon 
spectrum obtained using propagators smeared at both source and sink, at only a small computational cost. We 
investigate the matrix of correlators constructed from local and smeared operators, and are able to expose excited 
states of both the mesons and baryons. 



1. Introduction 

In this talk we present the results of an in- 
vestigation into gauge invariant smearing on a 
24'^ X 48 lattice at /3 = 6.2 with a hopping param- 
eter K = 0.152. This choice of k corresponds to 
a pseudoscalar meson mass of approximately 600 
MeV if we take the scale from the string tension 
(a-i ~ 2.73(5) GeV §). We begin by compar- 
ing effective masses obtained from local correla- 
tors with those from Wuppertal or "scalar wave 
function" smeared correlators . We introduce a 
variant of the iterative Wuppertal scheme which 
enables us to achieve the results of scalar wave 
function smearing but at a much lower cost in 
computer time. Using this scheme, we perform 
detailed measurements of the hadron spectrum 
using local and smeared sources and sinks. Fi- 
nally, we investigate the 2x2 matrices of cor- 
relators formed from propagators computed with 
these four source-sink combinations, and by di- 
agonalising these matrices are able to isolate the 
first excited states. 
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2. Gauge Invariant Smearing 

2.1. Wuppertal Smearing 

Wuppertal smearing uses a scalar propagator 
to smear a delta function source at the origin . 
The smeared source is given by the solution of the 
three dimensional, gauge invariant, Klein- Gordon 
equation. For our initial investigation of the de- 
pendence of the effective masses on the smear- 
ing radius, quark propagators were computed on 
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Figure 1. A comparison between the effective 
mass of the nucleon obtained using LL and LS 
propagators. 



an ensemble of 7 configurations using Wuppertal 
smeared sources at ks = 0.180 and 0.184, cor- 
responding to r ~ 2 and r ~ 4 respectively 
The effective mass of the nucleon obtained us- 
ing LL and LS propagators on the same ensemble 
of configurations is shown in figure |l| (we adopt 
the convention of LS representing a propagator 
smeared at the source and local at the sink). The 
lightest state of the nucleon is isolated nearer the 
origin with increasing smearing radius. 

The Wuppertal smearing algorithm creates sig- 
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nificant computational overheads when used to 
smear the propagator at the sink; the smearing 
algorithm must be implemented on every time 
slice, not just the source time slice, and for ev- 
ery spin component. Indeed, for ks = 0.184, the 
computational effort required to smear at the sink 
is comparable with that required to compute the 
propagator, and larger values of the smearing ra- 
dius are prohibitively expensive. 

2.2. Jacobi Smearing 

An alternative smeared source, J{x), can be 
obtained by solving the three dimensional Klein- 
Gordon equation as a power series in ks, stopping 
at some finite power N. 

N 

J{x) = ^ K5A 5xfi (1) 

ri=0 

where 

3 

^1=1 

This can be achieved easily using Jacobi iteration. 
When Ks is smaller than some critical value, the 
power series converges and one recovers the scalar 
wave function. For sufficiently large Kg, the se- 
ries diverges but provides an acceptable smeared 
source for suitable choices of N. A similar it- 
erative scheme has been used by the Wuppertal 
group !§. 

We found that for a radius of approximately 4, 
smearing with the Jacobi iteration was a factor 
of ten faster. The remainder of this paper will, 
therefore, concentrate on results obtained using 
Ks = 0.250, N — 50, corresponding to r ~ 4, on 
the same 18 configurations which were used in |^]. 

3. Hadron Masses 

We present in table 1 the masses of the pseu- 
doscalar, vector, nucleon and A using the LL, LS, 
SL and SS propagators, together with the corre- 
sponding time ranges used for the fits. All fits are 
performed using the full covariance matrix, with 
the errors extracted using a bootstrap analysis 
Q . The improvement in the determination of the 



Table 1 

The estimates of baryon and meson masses 

LL(12-16) LS(11-16) SL(11-16) SS(7-12) 
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baryon masses through the use of smeared sources 
and sinks is substantial. We expect that the study 
of a larger sample of configurations would lead 
to a still more marked improvement through the 
ability to perform a fit to the full covariance ma- 
trix for the SS correlators over a larger fitting 
range than that quoted in table 1. The lack of 
improvement in the determination of the meson 
masses suggests that mesons are already well ex- 
posed using local sources and sinks. 

A recent analysis Q of the hadron spectrum 
obtained using both wall sources and Wuppertal 
sources suggests that there is a systematic differ- 
ence in the baryon masses in the two cases, and 
that this difference is particularly noticeable in 
the determination of the A mass. To test for such 
a discrepancy in our data, we perform a bootstrap 
analysis of the A mass differences for the various 
combinations of sources and sinks, using the time 
ranges of table 1. For the LS and SL correlators, 
we find 

- = 0.026 ±11 . (2) 

This discrepancy is a Icr effect. However, since 
the expectation values of the LS and SL correla- 
tors should be identical, we attribute the discrep- 
ancy to limited statistics. 

4. Matrix Correlators 

By studying the matrix of correlators formed 
using the local and smeared operators, Ol and 
Os respectively, we can attempt to extract 
masses for the first excited state. We define the 
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2x2 matrix, C{t), of timeslice correlators by 

"^W-y ^^^^(i) uj\ss{t) J 
with 

The factor uj arises from the arbitrary normalisa- 
tion of the smearing function. Inserting a com- 
plete set of states in equation (|^), we obtain 

oo — -E„t 

cvit) = 5](0|O.(0,0)|n)(n|O](0,0)|0)^^. (4) 

n— ^ 

We will now discuss two methods for studying the 
matrix C(t) §. 

Diagonalisation of the Transfer Matrix 

Consider the eigenvalue equation 

C{t)u^ X{t,tQ)C{to)u. (5) 

for fixed to- If the system comprises only two 
independent states, then the eigenvalues of eq. |^ 
are 

X+{t,to) = e-(*-*°)^° 

X-{t,to) = e-(*-*«)^\ (6) 

The two states are separated exactly, and the co- 
efficients in eq. ^ grow exponentially with to- In 
general, where there are more than two states, 
for sufficiently large t we expect two states to be 
dominant. The coefficients of the contributions 
of the higher states to eq. B| do not exhibit ex- 
ponential growth with to [|||. Hence ideally we 
wish to study X-{t,to) for to as large as possible. 
However, the increase in the noise in the data far 
from the source generally requires that we choose 
to close to the origin. 

Diagonalisation of C{t) 

We can compute the eigenvalues, x+{tj'^) and 
X-{t,uj), of C{t) directly, and at large times, t, 
obtain 

X+it,u;) = /+He-^«*{l+5+H0(e-^^^*)} 

X-{t,^) = /-Me^^^*{l + .9-M0(e-^^*)}, 

(7) 



where 

A_£/i = El — Eq 

AE = min[£;i -Eo,E2- Ei]. (8) 

As in the previous method, the eigenvalues 
X+{t,oj) and X-C^jI^) are dominated by the 
ground state and first excited state respectively. 
However, even when there are only two states, 
corrections to the leading behaviour of the eigen- 
values remain. The coefficients of these correc- 
tions depend on the arbitrary parameter uj and, 
since we diagonalise C{t) at each timeslice, these 
coefficients are not positive definite. We will ex- 
ploit the dependence on w to seek cancellations 
between contributions from higher states, and 
hence extend the plateau region in the effective 
mass closer to the source. 

Such an approach has clear dangers. In par- 
ticular, the effective masses can approach their 
asymptotic values either from above or from be- 
low as we vary uj. One restriction does temper 
the uncontrolled nature of this method. We ob- 
serve that for sufficiently large times the excited 
state contributions to eq. |^ will be negligible, and 
thus the effective masses derived from x+i'^i'^) 
and x_(i,Ci;) must be insensitive to u. 

4.1. Results 

In figure ^, we show the effective mass derived 
from x~{^j ^) for the first excited state of the nu- 
cleon at several values of uj. We observe there is 
a region in t for which the effective masses coin- 
cide for all ijj. As noted above this can be taken 
as a signal that the first excited state has been 
isolated. The plateau region in the effective mass 
appears to extend closer to the source as uj de- 
creases. We attribute this to cancellations be- 
tween contributions of the higher excited states 
to x-(i,w). 

In figure ||we show the effective mass of the first 
excited state of the nucleon derived from A_ (t, to) 
with to = 1; and derived from X-{^j ^) at the op- 
timal value of UJ. The two methods yield consis- 
tent plateau, and using Method 2 we obtain 

mAr(4-8) = 0.98^ I , (9) 

where the fitting range is shown in brackets. 
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Figure 2. The variation with w of the effective 
mass derived from x_(t,a-') for the first excited 
state of the nucleon. AU values of lo are quoted 
relative to Wq, where Wq is defined as wf|j at 
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Figure 3. The effective mass of the first excited 
state of the nucleon derived from A_(t, to), and 
X-(^,^)- 



Comparing the computed ratio of the mass of 
the first excited state and that of the ground state 
for the nucleon with its physical value for light 
quarks, we find it is appreciably larger. We find 
mAr/mjv — 2.0, whereas TO7v(i44o)/™Af = 1-5- If 
this is disappointing, it should be noted that a 
mass of 0{a~^) is likely to be subject to con- 
siderable uncertainties arising from the non-zero 
lattice spacing. However, an unexpectedly large 
value for the mass of the nucleon excited state 
has been observed in other simulations 0. 
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